Nonlocal quark model beyond mean field and QCD phase transition 
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A nonlocal chiral quark model is consistently extended beyond mean field using a strict 1/N C expansion scheme. 
The parameters of the nonlocal model are refitted to the physical values of the pion mass and the weak pion decay 
constant. The size of the 1/N C correction to the quark condensate is carefully studied in the nonlocal and the 
usual local Nambu-Jona-Lasinio models. It is found that even the sign of the corrections can be different. This 
can be attributed to the mesonic cut-off of the local model. It is also found that the 1/N C corrections lead to a 
lowering of the temperature of the chiral phase transition in comparison with the mean-field result. On the other 
hand, near the phase transition the 1/N C expansion breaks down and a non-perturbative scheme for the inclusion 
of mesonic correlations is needed in order to describe the phase transition point. 



1. Introduction 

Understanding the QCD phase diagram is 
one of challenging issues in modern theoretical 
physics. A description of the most interesting re- 
gion of phase diagram at low and moderate tem- 
peratures/densities requires a non-perturbative 
approach, which also provides a proper under- 
standing of the chiral quark dynamics and the 
confinement mechanism. 

Until now, the only method which is directly 
based on QCD and which meets these require- 
ments is lattice gauge theory. Unfortunately, the 
application of lattice results to experimental data 
is complicated by the fact that most lattice cal- 
culations are performed with rather large quark 
masses, leading to unphysically large pion masses. 

In the present contribution we want to dis- 
cuss an effective model of low-energy QCD, ca- 
pable of describing the chiral as well as the de- 
confinement transitions. As a basis we use the 
PNJL model, see e.g. [1I2I3I4I5] . which general- 
izes the well-known Nambu-Jona-Lasinio (NJL) 
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model for the chiral quark dynamics by coupling 
it to the Polyakov loop, which serves as an order 
parameter of the deconfinement transition. 

To get a consistent picture of the hadronic 
phase it is important to go beyond the mean- 
field approximation and to include mesonic cor- 
relations. In the present work we suggest an im- 
provement of the PNJL model within strict 1 /N c 
expansion scheme and restrict ourselves to the 
case of zero chemical potential. 

2. Mean field 

The quark sector of the nonlocal chiral quark 
model is described by the Lagrangian 

C q = q(x)(iJp- m c )q{x) + ~[J${x) + J*(x)} , (1) 

where m c is the current quark mass, and = 
— iA^ the covariant derivative with a back- 
ground gluon field A M = A a ^ = S^qAq. The 
nonlocal quark currents are 

Jm(x) = Jd 4 (xix 2 )f(xi)f(x 2 )x 

xq(x - xi)T M q(x + x 2 ), (2) 
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where T a = 1, T n = ij 5 r a with a = 1, 2, 3. Spon- 
taneous breaking of chiral symmetry leads to the 
formation of a quark condensate and generates a 
dynamical contribution to the quark mass. As a 
result the Euclidean quark propagator takes the 
form 

S p = (i$ + m(p 2 )) -1 , m(p 2 ) = m c + m d f 2 (p 2 ) 

where f 2 (p 2 ) = exp(— p 2 /A 2 ) is a (Fourier trans- 
formed) Gaussian form factor and nid is an order 
parameter for dynamical chiral symmetry break- 
ing. The chiral condensate is obtained from the 
non-perturbative part of the quark propagator, 
57 = S p - S|, S c p = (ij) + m c )-\ 

External vector and axial vector fields can be 
introduced by a derealization of quark fields with 
the help of a Schwinger phase factor. Practically, 
nonlocal vertices with external currents can be 
obtained by using rules for derivatives of contour 
integrals. 

3. Beyond mean field 

Corrections to the dynamical quark mass be- 
yond the mean-field can be accounted for in a 
systematic 1/N C expansion scheme |6I7I8] of the 
quark selfenergy, = ij>A v + B p (see Fig. [TJ), 
and the quark propagator 

I ctmf+Nc\ _1 _ o-l 1 v Nc 
\O p ) — D p -f Z^ p 

crnf+Nc ^ q q Q _1_ /'Q't 

Jp ~ Op — OpL-ip Op -f ... \o) 

In order to arrive at a consistent approximation, 
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Figure 1 . 1 /N c corrections to the quark propaga- 
tor 



one needs to take into account 1/N C corrections 



to the meson propagator, see Ref. [5] for the NJL 
model and [5] for its nonlocal generalization. For 
the present model, we employ a diagrammatic 
technique developed in [TU]. The 1/N C correc- 
tions to meson properties will affect the results 
for the quark condensate via the readjustment of 
the model parameters (A, m c , GA 2 ) which are to 
be chosen such that the physical values for pion 
mass M^± — 139.57 MeV and weak pion decay 
constant f„ = 92.42 MeV are obtained at T = 0, 
while the dimensionless coupling GA 2 is left as 
a free parameter. Different parameterizations of 
the nonlocal model beyond mean field are given 
in Tab.[TJ The corresponding quark condensate is 
presented in Fig. [51 We will use the parametriza- 
tion No. 4 for finite T calculations because it pro- 
vides the highest (pseudo)threshold value of ex- 
ternal momentum in quark loop. 



Table 1 



Different parameterizations fitted to M^± and /, 



No. 


A, MeV 


m c , MeV 


m d , MeV 


GA 2 


1 


1479.2 


2.82 


139.2 


13.35 


2 


934.8 


5.58 


211.2 


14.89 


3 


705.9 


8.64 


269.1 


17.06 


4 


670.3 


9.31 


281.9 


17.64 


5 


580.5 


11.78 


322.5 


19.72 


6 


500.8 


14.95 


373.8 


22.83 


7 


445.3 


18.15 


424.0 


26.33 


8 


404.4 


21.37 


473.4 


30.20 



In table [2] we present the mean-field contribu- 
tions to pion mass and weak decay constant for 
different parameter sets together with values for 
fir estimated from Goldberger-Treiman and Gell- 
Mann-Oakes-Renner relations. For lower values 
of the current (and dynamical) quark masses the 
1 /N c corrections to pion mass and weak pion de- 
cay constant amount to 15 MeV and 20 MeV, 
respectively. For set 4 the corrections are only 
about 2 MeV and 5 MeV. 

Fig.[2]shows that the 1/N C correction to the ab- 
solute value of quark condensate is positive for all 
sets of model parameters. In the local NJL model 
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Figure 2. Quark condensate in the nonlocal PNJL 
model in MF approximation (dotted line) and in 
the 1/N C approximation with parameter adjust- 
ment (solid line) . The hatched region corresponds 
to the QCD sum rule limits for the quark conden- 
sate: 200 < -{qqy/ 3 [MeV} < 260. The vertical 
line corresponds to the dimcnsionlcss coupling for 
the parameter set No. 4 used in this work. 



of Ref. [8] it was found that this correction is neg- 
ative. However, in the local NJL model due to 
its non-renormalizability it is necessary to intro- 
duce different regularizations for the pure quark 
and the meson-quark loops, respectively. In [5] 
a Pauli-Villars regularization has been used for 
quark loops and a three-dimensional momentum 
cutoff Am for meson-quark loops. In order to 
study the transition from the local model to a 
local one let us construct a nonlocal model with 
three parameters 

1 . parameter of nonlocality A 

2. parameter of quark loop regularization A 9 

3. parameter of meson loop regularization Am 

The local model corresponds to the limit 

A^oo, A q = K p q hys , A M = A^ S , (4) 

while the nonlocal model without regularization 
can be obtained by setting 

A = A^ S , A 9 ^oo, A M ^oo. (5) 



Table 2 

Mean- field contributions to M n , / ff and estimates 
for f n from low-energy theorems. 



No. 




J 7T 


J 7T 


j-UMOH 

J 7T 


1 


155.5 


72.6 


93.6 


94.7 


2 


144.6 


83.4 


92.8 


91.3 


3 


142.5 


87.1 


92.0 


89.9 


4 


142.2 


87.6 


91.8 


89.6 


5 


141.7 


88.7 


91.4 


88.7 


6 


141.4 


89.6 


90.9 


87.6 


7 


141.2 


90.0 


90.6 


86.4 


8 


141.1 


90.3 


90.4 


85.3 



For definiteness, let us compare the local model 
[5] with the nonlocal from [5] with parametriza- 
tions fixed in the MF approximation. Note that 
for the given parametrizations, the MF quark 
condensates in the local and the nonlocal model 
agree within less than 0.5 %. 

The next step is to consider the 1/N C correc- 
tions and to investigate the role of the mesonic 3D 
cut-off Am- For this purpose it is very instructive 
to study the ratio of the full quark condensate to 
the MF contribution {qq)/{qq) mf . In Fig. [3] we 
compare the Am dependence of this ratio for the 
local NJL model as given in Ref. [5] (dotted line) 
to that of the nonlocal model (bold solid line) 
and its local limit (thin solid line). It is very in- 
teresting that in the region below 1.5 GeV these 
models predict a negative sign for the 1/N C cor- 
rection whereas for large mesonic cut-off the sign 
is positive. However, in the nonlocal model the 
absolute value of the correction saturates for Am 
larger than ~ 2.5 GeV, which is well above actual 
parametrizations for A g and Am in [Hj- 

4. Finite Temperature 

Details of the finite temperature extension of 
the nonlocal model can be found in [5110] . For the 
Polyakov loop potential we adopt the logarith- 
mic form of Ref. [4] which has been fitted to the 
quenched lattice data. In Fig. |4] we show the re- 
sulting temperature dependence of the quark con- 
densate (qq) T (normalized to its vacuum value) 
together with that of the Polyakov loop expec- 




Figure 3. The ratio (qq)/{qq) ml as a function of 
the meson cutoff Am- Local results (dotted line) 
is taken from Fig. 3a j8j. "Local" check (thin solid 
line) means the local limit of nonlocal calculations 
and thick solid line is nonlocal result. 



tation value $. In the nonlocal model without 
Polyakov loop the critical temperature for the chi- 
ral restoration is T c = 116 MeV, whereas the pure 
gauge sector has a critical temperature for decon- 
fincmcnt = 270 MeV, fixed from lattice data 
for $. When coupling the quark and gluon sec- 
tors, these critical temperatures get synchronized 
so that T c K,T d K, 200 MeV at the MF level. 

Near T c there is a wiggle in the behavior of the 
quark condensate which is caused by the pertur- 
bative nature of 1/N C expansion scheme. Namely, 
in the region of temperatures 183-223 MeV the 
correction to the quark condensate is larger than 
1/N C which one can naively expect. So, it seems 
reasonable to use some interpolation between the 
stable regions at low and high temperatures. 

On the other hand, 1/N C corrections slightly 
lower the temperature of the chiral phase tran- 
sition in comparison with the mean-field result. 
For low temperatures < 100 MeV our result for 
the quark condensate practically coincides with 
the xPT result, whereas the high T region is 
well controlled by the mean field. Near the phase 
transition the perturbative 1 /N c expansion breaks 
down and the prediction of the strict 1/N C ex- 
pansion scheme seems not reliable in a region of 
±20 MeV around the phase transition point. 



Figure 4. Temperature dependence of the quark 
condensate (thick solid line) and the Polyakov 
loop (thick dashed line) in the nonlocal PNJL 
model beyond mean field. Thick dash-dotted line: 
mean field contribution in the 1/N C model; thin 
dotted line: lowest order chiral perturbation the- 
ory (xPT); thin solid line: naive polynomial in- 
terpolation in the unstable region of the 1 /N c ex- 
pansion. 
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